Tiling the plane without supersymmetry 
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We present a way of tiling the plane with a regular hexagonal network of defects. The network is 
stable and follows in consequence of the three-junctions that appear in a model of two real scalar 
fields that presents Z3 symmetry. The Z3 symmetry is effective in both the vacuum and defect 
sectors, and no supersymmetry is required to build the network. 
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Domain walls appear in diverse branches of physics, 
envolving energy scales as different as the ones for in- 
stance in magnetic materials |Q] and in cosmology 
They live in three spatial dimensions as bidimensional 
objects that arise in systems with at least two isolated 
degenerate minima. In field theory they appear in the 
(3,1) dimensional space-time, and this may happen in 
supersymmetric theories, although supersymmetry plays 
no fundamental role for the presence of domain walls. 

Very recently, in a paper by Gibbons and Townsend 
1^, and also in Refs. [Q,^, one investigates the presence 
of domain walls and their possible intersections in a Wess- 
Zumino model, with a polynomial superpotential. In the 
supersymmetric theory, one can classify the classical solu- 
tions as BPS and non-BPS states, according to the work 
of Bogomol'nyi, and of Prasad and Sommerfield The 
BPS states are stable, and are expected to play some 
role in investigating duality in supersymmetric models. 
We recall that no BPS state can be annihilated under 
continuum variation of the parameters that define the 
supersymmetric theory. 

In Ref. one investigates models of coupled real scalar 
fields in bidimensional space-time. These investigations 
provide a concrete way of finding BPS states and suggest 
other studies, in particular on the subject of defects inside 
defects - see Ref. ||] . Most of the models investigated in 
|@,|| can be seen as real bosonic portions of supersym- 
metric theories. In supersymmetric models the presence 
of discrete symmetry may produce BPS and non-BPS 
defects. The BPS states lie in shorter multiplcts, and 
preserve the supersymmetry only partially [plJlCfl . There 
are BPS states that preserve 1/2 of the supersymmetry, 
but the possibility of BPS states preserving 1/4 super- 
symmetry is subtler, and is shown to appear as junctions 
|PT| , P^ of domain walls in the recent papers . 

In the present work we start dealing with the bosonic 
portions of supersymmetric theories. We do this guided 
by the discrete Z3 symmetry, with the aim of describing 
the presence of three-junctions and the network of defects 
that it can generate. We first point out that supersymme- 
try introduces restrictions that may lead to instability of 
the junction, or at least of the network that it could gen- 
erate. We then examine another model, and show that 
all the difficulties found in the supersymmetric context 
are circumvented by just giving up supersymmetry. 



The subject of this work may be of interest to several 
different branches of physics, in particular in applications 
concerning the entrapment of networks of defects inside 
domain walls. This possibility can be implemented with 
three scalar fields, in models engendering the Z2 x Z3 
symmetry, following the lines of Refs. Other ap- 

plications may include strong interactions, if we recall 
that the group is the center of the SU (3) group - see 
for instance Refs. |ljjl^. Also, there are applications 
to systems of condensed matter, in particular on issues 
concerning pattern formation ]l6t , as for instance in the 
case of the thermal convection studied in Ref. p7[ |. 

We start describing the two real scalar fields (/) and x 
in bidimensional space-time with the Lagrangian density 



V 



(1) 



Here V = F(0, x) is the potential. In the supersymmetric 
case it has the general form 
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where W — W{<l),x) is the superpotential. 

The superpotential allows introducing several proper- 
ties, as shown in Refs. 0,^. For instance, for static 
fields the equations of motion are solved by first-order 
differential equations d(j)/dx = and dx/dx = W^. 
The energy of solutions of the first-order equations are 
given by E% = lAI^y l, with Al^y = Wi - Wj and 
Wi — W{(j)i,Xi), where {(j)i,Xi) represents the i-th vac- 
uum state of the model. This is the Bogomol'nyi bound, 
and the corresponding solutions are BPS solutions. The 
BPS solutions are linearly or classically stable. 

We guide ourselves toward the topological solutions by 
introducing the topological current 



J" 



d0 



(3) 



It obeys 9qJ" = 0, and it is also a vector in the ((l),x) 
plane. For static configurations we have J" = p* p, 
where p = p((j), x) is the charge density. This charge den- 
sity allows writing p*p as twice the kinetic energy density 
of the topological solution, and this can be used to infer 
stability of junctions - see below. 
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Let us now consider a specific model, defined by the 
superpotential 




A X + 



The first-order differential equations are given by 
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The model is described by the fourth-order polynomial 
potential 



+3X 






+2A^ X 
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Thus, it behaves standardly in one, two, and three spatial 
dimensions. 

The above potential presents the three vacuum states 
(0, v/473), (-1, - 073): and (1, - v/l73)- These min- 
ima form an equilateral triangle invariant under the 
rotations, with side I — 2. The values of the super- 
potential at the minima are Wi — —A, W2 — 2X and 
W3 = -2A. The energies of the BPS states are |A|, 3 |A|, 
and 4 |A|. All the three sectors are BPS sectors, but they 
do not present the Z3 symmetry that connects the vac- 
uum states. The highest energy is associated to the sector 
connecting the second and third vacuum states. This is 
the only sector where we can find explicit solutions. They 
are given by (j}{x) = — tanh(3Aa:) and xi^) = — \/l/3 . 
This is a BPS state, representing an orbit in the {4>, x) 
plane. The orbit is a straight-line segment that connects 
the corresponding vacuum states. The orbits connecting 
the other vacua cannot be straight-line segments. They 
cannot be obtained by rotating the (0, x) plane accord- 
ing to the Z3 symmetry, and so the defect sectors do 
not present the Z3 symmetry that connects the vacuum 
states. This fact also appears when one identifies the 
tensions of the BPS defects. They are given by ti = |A|, 
t2 = 3 |A|, and ^3 = 4 |A|. They are different and do not 
obey the Z3 symmetry. They are such that — ti + t2, 
and do not strictly obey the triangle inequality one needs 
to ensure stability ||T^ of the three-junction that appears 
in this model. 



To circumvent instability of the three-junction we now 
follow Refs. We make contact with these works after 
considering superpotentials that satisfy W^^ + W^^ — 0. 
In this case, for harmonic superpotentials one adds to the 
two first-order eq uations d(j)/dx — and dx/dx — 
the two new [|l3| first-order equations: d(p/dx — —W^ 
and dx/dx = W^. Solutions to these equations also min- 
imize the energy and solve the equations of motion. This 
allows introducing W^(0, x) such that W(j, = —W^ and 
= W^. We use W and W to introduce the complex 

superpotential, W = W + iW. We write the complex su- 
perpotential in terms of the complex field (p+ix, and this 
is the way one gets from the investigations of Refs. to 
the recent possibility of describing three-junctions 

preserving 1/4 supersymmetry. However, junctions re- 
quire the presence of at least three minima, and this is 
only achieved when the superpotential is of at least the 
fourth-order power in the complex field. This means that 
the model behaves standardly only in one and two spa- 
tial dimensions. In this case one can show explicitly 
that the three-junction is stable and breaks 1/4 super- 
symmetry, although supersymmetry itself does not allow 
the presence of a stable network of defects ■ Owing 
to the fact that each adjacent junction in the network 
has opposite winding number, any adjacent vacua should 
be connected with defect solutions also having opposite 
winding numbers along the same orbit |^ . Since we have 
to use different conjugate Bogomol'nyi equations to take 
into account these winding numbers, the network cleary 
cannot be BPS and then can decay. 

We then give up supersymmetry, turning our attention 
to polynomial potentials that engenders the Z3 symme- 
try, and that supports stable three-junctions that gener- 
ate a regular hexagonal network of defects. Interestingly, 
we have found a fourth-order polynomial potential that 
do the job. It is given by 

mx) = AV ('^'-^)+AV {x'~l 

+ 2\'4>'x'~X'4>{^'~ix') + ^{x' (8) 

This potential was introduced in Ref . ||l8| . The equations 
of motion for static configurations are 
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dx 
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The potential has three degenerate minima, at the points 
vi = (3/2) (1, 0) and V2,3 = (3/4) (-1, ±V3). These min- 
ima form an equilateral triangle, invariant under the Z3 
symmetry. The distance between the minima is (3/2)-\/3. 

We can obtain the topological solutions explicitly. The 
easiest way to do this follows by first examining the sector 
that connects the vacua V2 and W3. This is so because in 
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this case we set (p = —3/4, searching for a strainght-Hne 
segment in the (0, x) plane. This is compatible with the 
Eq. and reduces the other Eq. ( p!o| ) to the form 



(11) 



This implies that the orbit connecting the vacua V2 and 
f3 is a straight line. It is such that, along the orbit the 
X field feels the potential [x^ - (27/16)]2. This shows 
that the model reduces to a model of a single field, and 
the solution satisfies the first-order equation 



dx 



= V2X 
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The solution is 



x(x) — — --\/3tanh 




(13) 



The other solutions can be obtained by rotations obeying 
the Z3 symmetry of the model. 

The full set of solutions of the equations of motion are 
collected below. In the sector connecting the minima V2 
and V3 they are 



^(2.3) 
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In the sector connecting the minima vi and V2 they are 

,(±1 ■J , y , , / /27 
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y.-i — — ± — tanh 
(1,2) 8 8 



(±) 3 
^(1,2) 8 ■ " ' 8 
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In the sector connecting the minima vi and V3 they are 

f27 
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^(S) = "^^"F ^%/3 tanh ^y^Aa;^ 



(18) 
(19) 



The label (±) is used to identify kink and antikink. All 
the solutions have the same energy, (9/4) -\/27/8 |A|. 

We examine how the bosonic fields behave in the back- 
ground of the classical solutions. We do this by consid- 
ering fluctuations around the static solutions (t>{x) and 
x{x)- We use the equations of motion to see that the 
fluctuations depend on the potential 



U(x) 



(20) 



Evidently, after obtaining the derivatives we substitute 
the fields by their classical static values (j){x) and x{^)- 
The model under consideration is defined by the potential 
(^). In this case we use ( |l^ and (^5|) to obtain two 
decoupled equations for the fiuctuations. The potentials 
of the corresponding Schrodinger-like equations are 



27 

Unix) = — A^ 
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27, 



2 sech^ 



4 — 6 sech 




(21) 



(22) 



(12) The eigenvalues can be obtained explicitly: in the x di- 
rection we get Wq=0 and — (9/2)-y/A^/2, and in 
the direction we 

have = (9/2)v/AV2- This shows 
that the pair ( p^ ) and (^5|) is stable, and by symmetry 
we get that all the three topological solutions are stable 
solutions. 

The classical solutions present the nice property of hav- 
ing energy evenly distributed in their kinetic (k) and po- 
tential (p) portions. In terms of energy density they are 



k{x) 



( ^ 1/^27 
P^") = 4 ¥ 



A^ sech* 




(23) 



To understand this feature we recall the calculation done 
explicitly in the sector with = —3/4, constant. There 
the model is shown to reduce to a model of a single field, 
a model that supports BPS solutions. Within this con- 
text, the above solutions are very much like the non-BPS 
solutions that appear in supersymmetric systems Jll3[ . 
We use this property and the topological current (Ifto 
obtain p*p = e, where e{x) = k{x) + p{x) is the (total) 
energy density of the solution. We use this result and the 
notation ij, to identify the sector connecting the vacua 
{4>i,Xi) E^nd {(j)j,Xj): to show that for any two different 
sectors ij and jfc, i, j,k — 1,2,3 we get that 



{Ptj + Pjkf {Pij + Pjk) < pljPtj + p]kP]k 



(24) 



This condition shows that the three-junction is a process 
of fusion of defects that occurs exothermically, providing 
stability of junctions in the present model. This result 
is more general than the one in Ref. |p^ , which appears 
within the context of supersymmetry. Evidently, our re- 
sult also works for BPS and non-BPS solutions that ap- 
pears in supersymmetric models, with the property of 
having energy evenly distributed in their kinetic and po- 
tential portions ||l^ . 

We notice that the orbits corresponding to the stable 
defect solutions form an equilateral triangle in the ((/), x) 
plane. This is so because the solutions are straight-line 
segments joining the three vacuum states in configuration 
space. They are degenerate in energy, and this allows 
associating to each defect the same tension 



9 /27 
t = - \ — A 
4 V 8 ' ' 



(25) 
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This makes < tjk +tki, i, J, fc = 1, 2, 3, and now the in- 
equahty is strictly vahd in this case, stabiUzing the three- 
junction that appears in this model when one enlarges the 
space-time to three spatial dimensions. 

We consider the possibility of junctions in the plane, 
which may give rise to a planar network of defects. We 
work in the (2,1) space-time, in the plane {x,y). We 
identify the plane {x, y) with the space of configurations, 
the plane (0, x)- We illustrate this situation by consider- 
ing, for instance, the solutions we have already obtained. 
They are collected in Eqs. (P^)-(]l9|) in (1,1) dimensions. 
In the planar case they change to 



»(2,3) 



3 

'4 




and 



^1,2) 



= ±^%/3tanh 



3 9 , I 27 , , /- ' 
- T g tanh I 2 Y ''^ ^ ' 



= ± ^%/3tanh (^l^m^iy + V3x) j 



and 



^(1,3) 



3 9 , / 1 /27 , , /- ' 
- ± - tanh I -y — A (y - V3x) 



(26) 
(27) 

(28) 
(29) 

(30) 



X(S)=-^^±^^tanh(^iy|A(2/-y3x)^ (31) 

These planar defects are domain walls, and can be used 
to represent the three-junction in the limit of thin walls. 




FIG. 1. A regular hexagonal network of defects, formed 
by three-junctions surrounded by domains representing the 
vacua 111 — 1,V2 = 2, and 113 = 3. 



The three-junction that appears in this Za-symmetric 
model allows building a network of defects, precisely in 
the form of a regular hexagonal network, as depicted in 
FIG. 1 in the thin wall approximation. In this network 
the tension associated to the defect is the tipical value 
of the energy in this tiling of the plane with a regular 
hexagonal network, which seems to be the most efficient 
way of tiling the plane. As we have shown, our model be- 
haves standardly in (3, 1) dimensions. It supports stable 
three-junctions that generate a stable regular hexagonal 
network of defects. 
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